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ABSTRACT 

We have compared the optical emission line ratios of type 2 quasars from Zakamska et al. 
with standard AGN photoionization model predictions, Seyfert 2s, HII galaxies, and narrow 
line FRII radio galaxies. Moderate to high ionization narrow line radio galaxies and Seyfert 2s 
are indistinguishable from type 2 quasars based on their optical line ratios. The standard AGN 
photoionization models, widely discussed for other type 2 AGNs, can reproduce successfully 
the loci and trends of type 2 quasars in some of the main diagnostic diagrams. These models 
are not exempt of problems and the discrepancies with the data are the same encountered for 
other type 2 AGNs. As for these, realistic models must take into account a range of cloud 
properties, as widely demonstrated in the literature. 

The Zakamska et al. sample is strongly biased towards objects with high line luminosities 
(L[OIII]>10 42 erg s" 1 ). We have found that stellar photoionization is obvious in a fraction of 
objects (3 out of 50) which are characterized by low [OIII] luminosities compared with most 
type 2 quasars in the sample. We suggest that if the sample were expanded towards lower 
line luminosities (^10 40 ~ 42 erg s _1 ) stellar photoionization might be evident much more 
frequently. 

We explore an alternative scenario to pure AGN photoionization in which a varying con- 
tribution of stellar ionized gas contributes to the line fluxes. Although the models presented 
here are rather simplistic and not strong quantitative results can be extracted regarding the rel- 
ative contribution of stellar vs. AGN photoionization, our results suggest that adding a varying 
contribution of stellar photoionized gas works in the right direction to solve most of the prob- 
lems affecting the standard AGN photoionization models. The "temperature problem", on the 
other hand remains. 

Key words: galaxies: active; quasars: general; quasars:emission lines 



1 INTRODUCTION 



Type 2 active galactic nuclei (AGNs) are those AGNs whose per- 
mitted and forbidden lines have similar values of full width at half 
maximum (FWHM). In many cases (e.g. Seyfert 2 galaxies) this 
can be explained as a consequence of obscuration of the central 
region by large column densities of gas and dust. In the standard 
unification model the obscuring structure is toroidal in shape so 
that the view to th e inner nuclear region is blocked for some orien- 
tations ( Antonucci 1 19931) . According to this model, certain classes 
of type 1 and type 2 AGNs are the same entities, but have differ- 
ent orientations relative to the observer line of sight. If this model 



is valid for the most luminous AGNs (quasars), there must exist a 
high-luminosity family of type 2 quasars. 

The existence of such an object class was predicted a long 
time ago, but it has been only in the last few years that type 2 
quasars have been discov ered in large quantitie s in X -ray, mid-IR 
(Martmez-Sansigre et al. 120051 . Szokoly et al. |2004) a nd optical 
surveys (Zakamska et al. 2003). Zakamsk a et al. I J2003I) identified 
~145 objects in the redshift range 0.3< z <0.8 in the Sloan Dig- 
ital Sky Survey (SDSS; York et al. l2000h with the high ionization 
narrow emission line spectra characteristic of type 2 AGNs and nar- 
row l ine luminosities typical of type 1 quasars (see also Reyes et al. 
2008). Their IR and X-ray properties are consistent with their inter- 
pretation as powerful obscured AGN (Ptak et al. 2006, Zakamska 
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et al. l2004b . They show a wide range of X-ray luminosities and ob- 
scuring column densities. About 40 objects in their sample were 
detected with IRAS and have infrared luminosities among the most 
luminous quasars at similar redshift. The host galaxies are ellipti- 
cals, although with irregular morphologies, and the nuclear optical 
emission is highly polarized (Zakams ka et a l. uOOq) . The detecti on 
rate in radio (~10%, Zakamska et al. l2004 Vir Lai & Ho l2007u is 
consistent with that of other AGN types. 

Spectroscopic studies of type 2 quasars have focused so far on 
identifying emission lines, measuring some basic parameters (line 
luminosities, redshift, line widths) and searching for correlations 
among them and with other observables (equivalent widths, color 
magnitudes, radio luminosities, etc). All this is critical to classify 
the objects, to investigate the nature of the powering mechanism, 
the obscuring structure and, ulti mately , test the validity of the uni- 
fication models (e.g. Reyes et al. l200a) . 

However, little work has been done to characterize the gaseous 
and ionization properties of type 2 quasars. Given this lack of 
knowledge, the primary goal of the work presented here is to use 
the emission line information to study the excitation mechanism, 
the physical conditions and ionization properties of the gas. 

Detailed spectroscopic studies of other active galaxy types, 
such as narrow line radio galaxies (NLRGs, e.g. Robinson et al. 
Il987l) have allowed during the last few decades the characterization 
of the chemical abundances of the gas, the ionization mechanism, 
the physical and kinematic properties, etc. Similar work must be 
done for type 2 quasars. These studies will ultimately provide valu- 
able information about the formation process of the host galaxy, 
the star forming and chemical enrichment histories and t he ori- 
gin of the nuclear activit y (e.g. Tadhunter , Fosb ury & Quinn 1 19891. 
Villar-Martm et al. l2005t Humphrey et al. l2008l) . An important dis- 
advantage of NLRG over type 2 quasar studies is that the radio 
activity, via jet-gas interactions, imprints important distortions on 
the properties of the ionized gas making it difficult to investigate 
the intrinsic properties of the host galaxy and environment, as well 
as the chemical abundances and the nature of the excitation mecha- 
nism (e.g. Tadhunter l 20020 . The study of radio-quiet type 2 quasars 
should not suffer from such effects. 

Throughout this paper we assume Q,a = 0.73, fl m = 0.27 and 
Ho = 71 kms" 1 Mpc -1 . 



2 THE OBJECT SAMPLE 

The objects studied in this paper are a sub-sample of candidate type 
2 quasars from t he Slo an Digital Sky Survey (SDSS) selected by 
IZakamska etHI d2003l) in the redshift ran ee 0.3< z <0.8 (~145 
objects). The selection criteria applied by Zakamska et alj d2003h 
are listed below (notice that not necessarily all criteria apply to all 
objects. See Zakamska et al. 2003 for more detailed information))^ 

• i < 19.1 

• 0.3 < z < 0.83 

• Signal to noise ratio >7.5 

• Equivalent width (EW) of [OIII] > 4 A 

• Luminosity of [OIII], L[OIII]>3 x 10 s L 

• Full width at half maximum, FWHM(H/3)<2000 km s _1 



1 Throughout the paper the emission lines will be named as follows: [OIII] 
for [Oin]A5007; [On] for [OII]A3727; [Nil] for [NII]A6583; [NeV] for 
[NeV]A3426; [Nelll] for [NeIII]A3689; Hell for HeIIA4686. 



• Selection of active galaxies among emission line galaxies, 
based on line ratio criteria (different line ratios depending on z) 

• Other signs of AGN activity such as the detection of the high 
ionization lines [NeV]AA346,3426 and/or FWHM([OIII])>400 km 

s" 1 . 

The sub-sample studied here contains 50 type 2 quasars, 
which have been selected to span the full z range and the full [OIII] 
luminosity range (^3x 10 s Lq) of the original type 2 quasar sam- 
ple. We did not set constraints on the line equivalent widths, al- 
though the selection criteria on the original sample do contain an 
EW criterion, as stated above. 

The spectra were corrected for Galactic extinction. A galaxy 
template spectrum was subtracted for objects with low line EWs, 
to correct for possible underlying stellar and interstellar absorption 
(Zhang, Dultzin-Hacyan & Wang 120071) . This procedure was nec- 
essary for a small fraction of objects (~15% only). 

The spectra were not corrected for internal dust reddening be- 
cause such correction was not possible for all objects. It should not 
affect our conclusions since we will compare our results with pre- 
vious works on other AGN types (Seyfert 2s, radio galaxies), in 
which no internal extinction correction was applied either. 



3 ANALYSIS AND RESULTS 

We investigate in this section the dominant ionizing mechanism of 
the optical line emitting gas in the type 2 quasar sub-sample: AGN 
vs. stellar photoionization. 

We will ignore shocks in our discussion a s an al ternative ion- 
ization mechanism (e.g. Dopita & Sutherland 1 19961) . In our sub- 
sample, 44 out of the 47 objects for which radio information is 
available are radio quiet (i.e. Liaghz <10 31 erg s _1 Hz -1 sr _1 ) 
and therefore, shocks induced by the radio structures are not ex- 
pected to play a significant role in t he ion ization of the gas (e.g. 
Clark et al. 1998, Villar-Martm et al. ll999h . In the vast majority of 
the type 2 quasars in the sample considered here, the bulk of the 
line profiles is characterized by rather quiescent kinematics, rather 
than perturbed, as one would expect if shocks were present. 

3.1 AGN photoionization 

We show in Fig. 1 several diagnostic diagrams involving optical 
emission lines in which we plot the location of the SDSS type 2 
quasar sub-sample (green, solid circles and blue solid triangles). 
For comparison, we plot in the same diagrams the lo cus of HII 
galaxies from the catalogue of iTerlevich et alj dl99 ll) (magenta 
small symbols). These are absent in diagrams involving the Hell 
and [NeV] because such lines are rarely detected in this object 
class. 

The solid lines represent the standard sequence of pho- 
toionization models, built wi th the multipurpose code MAPPINGS 
Ic (Binette, Dopita & Tuohv ll985l : Ferruit et al. ll997l) that repro- 
duces some of the main properties of the emission line spectra 
of narrow line radio galaxies ^it different redshifts (e.g. Robinson 
et al. Il987l Humphrey et al. 120081) . The ionizing continuum is a 
power law of index a=1.5 {F u oc v~ a ), with a cut off energy of 

2 U = j - J2 n — -, where Q is the ionizing photon luminosity of the 
source, r is the distance between the cloud and the ionizing source, nu is 
the hydrogen density at the illuminated face of the cloud and c is the speed 
of light. 
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Power law (alpha —1.5) 

• HII galaxies (Terlevich et al. 1991) 

• Type 2 quasars 

▲ Hybrid Type 2 quasars 



Figure 1. Loci of type 2 quasars in diagnostic diagrams involving the main optical emission lines detected in most spectra. Hybrid objects are those with 
signatures of stellar photoionization. The solid line is the standard sequence of solar metallicity AGN photoionization models along which the ionization 
parameter U varies. \g(U) values are shown for some models in several diagrams. The arrows indicate the sense of increasing U . Errorbars are shown in those 
diagrams where the measurement errors could have an impact on our interpretation and conclusions (d and g). In all other diagrams the errorbars are in general 
smaller. The comparison with the models and HII galaxies suggests that in general, AGN photoionization plays an important role in ionizing the gas and the 
ionization parameter varies from object to object. On the other hand, the standard AGN sequence shows some discrepancies with the data (e.g. diagrams d, e, 
/, g) which are the same encountered for other type 2 AGNs. Stellar photoionization is apparent in a small fraction of objects (hybrid object). 



50 keV. The clouds are considered to be isobaric, plane-parallel, 
dust-free ionization-bounded slabs of density n =100 cm at 
the illuminated face and characterized by solar abundances (An- 
ders & GrevesseJ989). Since the SDSS spectra integrate the emis- 
sion from nuclear (i.e, the narrow line region, NLR) and (possibly) 
extended gas, a range of cloud densities is expected with values 
as high as 10 6 cm -3 or more. Given the low critical density of 
[OII]A3727 (~3000 cm" 3 ), such high densities will pose difficul- 
ties to reproduce the observed strength of [Oil] relative to other 
emission lines, which suggest a substantial contribution to the in- 



tegrated spectrum of a low density component. Thus, although it 
is clear that assuming a single density is a rough simplification 
(characteristic on the other hand of the standard AGN sequence 
discussed in numerous papers), we will assume n e =100 cm -3 as 
a representative density of the gas within the SDSS aperture and 
mention the possible impact of higher densities when necessary. 

There are several diagrams which are specially sensitive to the 
ionization parameter U, in the sense that for a specific ionization 
mechanism, the line ratios involved are strongly influenced by a 
change in U. These diagrams are a, b, c, e and / (Fig. 1). The mod- 
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Figure 2. Comparison among type 2 quasars, Seyfert 2s and narrow line 
FRII radio galaxies (nuclear and extended emission line regions, EELR). 
The solid line is the AGN model sequence, as in Fig. 1. The three object 
classes overlap in the area of the diagram with lg([/)^-2.6. The lack of 
type 2 quasars with lg([OIII]/H/3) <0.5 is due to the selection criteria. 
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Figure 3. Type 2 quasars and HII galaxies plotted in Baldwin, Philips & 
Terlevich (1982) diagnostic diagram to classify emission line extragalactic 
objects. Type 2 quasars are clearly separated from HII galaxies. They oc- 
cupy the area of the diagrams where typical type 2 AGNs lie. Symbols and 
lines as in Fig. 1. 



els shows the large shift along the sequence expected as U varies. 
Looking at the distribution of data points in diagrams a, b, c and 
e, type 2 quasars follow a sequence whose shape can be accounted 
for by a variation in U. Most objects have lg(U) values in the range 
[-2.7,-1.7] 

There are several arguments that suggest that AGN photoion- 
ization plays an important role in ionizing the gas: 

• The shape and location of the sequence defined by type 2 
quasars is adequately accounted for by the AGN models in the a, b 
and c diagrams. 

• Type 2 quasars have similar line ratios than Seyfert 2 and nar- 
row line radio galaxies. We compare in Fig. 2 the loci of type 2 
quasars, Seyfert 2 (black open triangles) and radio galaxies (red di- 
amonds for the extended emission line regions (EELR) and red 'x' 
for the nuclear regions) in the [OIII]/H/3 vs. [OII]/[OIII] diagnostic 
diagram. This diagram is chosen because of its high sensitivity to 
the ionization level of the gas, whose variation, as proposed above, 
can explain the trends defined by type 2 quasars in those diagnos- 
tic diagrams with the highest sensitivity to U. The data for most 
Seyfert 2s were collected from lKoskilfl978h . Some more line ratios 
were taken from papers on individual sources (e.g. B ennert et al. 
I2006L Ferru it et al l 19991 , Contin i et al |2002|. Allen et al. ll999L Diaz 
etal.[l988, Kraeme r & Cr ensaw l200d) . The radio galaxy data set is 
the same as in lHoitl(f2006). None of these line ratios have been cor- 
rected for internal reddening. Fig. 2 shows that type 2 quasars over- 
lap in the diagram with radio galaxies with moderate and high ion- 
ization level (lg((7)>-2.6) and with most Seyfert 2 galaxies. This 
suggests similar ionization mechanism, ionizing continuum shape, 
physical properties of the gas and a similar range of chemical abun- 
dances. The lack of type 2 quasars with lg([OIII]/H/3) <0.5 is due 
to the selection criteria (see Zakamska et al. 2003). 

• Type 2 quasars are in general clearly separated from HII galax- 
ies (see diagrams a, b and d). This is particularly evident in dia- 
gram d, involving the [NII]/Hq ratio (see also Fig. 3). This is also 
evident in Fig. 3, which shows one of the most efficient diagrams 
at separating active galaxies and HII galaxies (Baldwin, Philips & 



Terlevich [l982l) . The [NII]/Ha ratio was measurable only for 30% 
of the objects. In a few cases, the lines are in a very noisy part of 
the spectrum due to strong sky residuals. For most spectra, they are 
outside the observed spectral range. In spite of this, the separation 
from HII galaxies (also with scarce [NII]/Hq measurements) is evi- 
dent and type 2 quasars occupy the area of the diagram where active 
galaxies lie (see Fig. 5 in Baldwin, Philips & Terlevich 1982). The 
harder continuum in active galaxies, produces a much deeper par- 
tially ionized region in which low ionization lines, such as [Nil] are 
efficiently excited. Although less evident, a similar effect accounts 
for the fact that the sequence defined by HII galaxies in diagrams a 
and b is shifted down relative to that defined by the type 2 quasars. 

• [NeV] and/or Hell are detected in all objects, while these lines 
are rarely found in the spectra of star forming galaxies. 

• The line luminosities of type 2 quasars are character istic o f 
quasars rather than of star forming objects (Zakamska et al. l2003l) . 

In spite of the general success of the standard AGN pho- 
toionization models to explain the emission line spectra of type 2 
quasars, the models encounter some problems, which are exactly 
those found i n studi es of other type 2 AGNs (e.g. Groves, Dopita 
& Sutherland |2004 Binette, Wilson & Storchi-Bergmann 1 1994 
Robinson et al. ll987l) . 

A. High ionization lines imply higher U values than low ion- 
ization lines 

The diagrams in Fig. 1 show that low and high ionization lines 
require different values of U. This is specially clear in diagrams e 
and / involving the [NeV] line, where type 2 quasars are shifted 
relative to the AGN model sequence. While the [NeV] ratios im- 
ply lg(f/) in the range [-1.8,-1], the implied range of values in dia- 
grams a, b and c is [-2.7,-1.6]. As for other AGN types at different 
redshifts ( e.g. H umphrey et al. [2008, Binette, Wilson & Storchi- 
Bergmann 1996), the models that reproduce some of the main line 
ratios under-predict the highest ionization lines ([NeV]A3426 in 
our case). 

Since the faint [NeV] line is detected in most objects plotted 
in Fig. 1, the discrepancy is not the result of different diagrams 
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Figure 4. The [OIII]AA5007,4959/[OIII]A4363 temperature sensitive ratio 
predicted by the AGN photoionization models is well above the measured 
values. Errors cannot account for such large discrepancy. Arrows indicate 
lower limits. Symbols and lines as in Fig. 1. 



showing different sub-samples. The problem is found in individual 
objects. 

This cannot be accounted for by reddening effects. Cor- 
recting for reddening would produce higher [NeV]/[OII] and 
[NeV]/[NeIII] that would result on larger U values. On the other 
hand, [OII]/[OIII] would also become larger, implying lower U val- 
ues. Therefore, the discrepancy would be worse. 

B. Too low electronic temperatures 

We have plotted in Fig. 4 the temperature sensitive 
[OIII]AA5007,4959/[OIIi]A4363 ratio vs. [OIII]/H/3. The arrows 
correspond to objects for which [OIII]A4363 was not detected and 
only upper limits could be estimated. For all objects the predicted 
and measured values of [OIII]AA5007,4959/[OIII]A4363 are dis- 
crepant by at least a factor of 2 (a shift of >0.3 in lg). This cannot 
be accounted for by measurement errors (see Fig. 4) or redden- 
ing effects, since correcting for this would result on even lower 
[OIII]AA5007,4959/[OIII]A4363 ratios. 

Higher densities can explain low 

[OIII]AA5007,4959/[OIII]A4363 values consistent with the 
data, but produce strong discrepancies with other line ra- 
tios. As an example, models with n=10° cm -3 produce 
lg([OIIl]AA5007,4959/[OIII]A4363K2, as measured for many 
type 2 quasars. However, for such models [OII]/[OIII]^0.02, 
[OII]/H/3 ^0.3 (much lower than the measured values, Fig. 1) 
and [OIII]/H/3 ^18 (higher than measured, Fig. 1). High densities 
only, therefore, in general do not solve the problem. 

In the low density limit (Osterbrock 1989), this discrepancy 
implies that the models predict too low electron temperatures: 
T e ^ 1 1000 K for all models of the U sequence, while the measured 
line ratios imply T e >15 000 K, being >20 000 in several cases, 
errors considered. The same problem has been discussed in detail 
for ot her type 2 AGNs (e.g. Binette, Wilson & Storchi-Bergmann 
1 19961 Robinson et al. ll987h . 

C. Too small scatter of the Hell/HP ratio 

The data present a large HeII/H/9 scatter inconsistent with the 
standard AGN sequence (Fig. 1 bottom panel), which is not due 
to errors in the measurements, neither reddening effects. There are 
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Figure 5. Examples of objects with very high (top) and very low (bottom) 
HeII/H/3 ratios ratios. 



objects with too high and objects with too low HeII/H/3 ratios com- 
pared with the model predictions (see Fig. 5). Such large scatter 
has been observed also in low z radio galaxies (e.g. Robinson et al. 
1987). 

D. Too strong [Nil] emission 

The [NII]/Hq vs. [OII]/[OIII] diagram (Fig. 1) shows a very 
large scatter in the [NII]/Ha ratio inconsistent with the standard 
AGN (solar metallicity) predictions which again cannot be ex- 
plained by errors in the measurements or reddening effects. For 
a large fraction of objects, the [Nil] emission is too strong. Given 
that this ratio is a direct metallicity indicator, rather than being a 
problem for photoionization models, the large range of values sug- 
gests that the nitrogen/hydrogen ratio varies substantially within 
the sample. Nitrogen is likely to be overabundant in those objects 
with large [NII]/Ha values (e.g. Robinson et al. 



3.2 Alternative AGN scenarios. 

Possible solutions to these problems have been extensively dis- 
cussed in the literature for more than 20 years. We present here 
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a brief summary of the main alternative AGN photoionization sce- 
narios. 

• Mixture of matter and ra diatio n bounded clouds (e.g. Binette, 
Wilson & Storchi-Bergmann 1 19961 : Viegas & Prieto ll992l) . Con- 
trary to ionization bounded (IB) clouds, matter bounded (MB) 
clouds are not sufficiently thi ck to absorb all the ionizing radiation. 
In the scenario proposed by iBinette. Wilson & Storchi-Bergmannl 
( 1996) the IB component sees a spectrum modified by absorption 
within the MB component. 

• Locally optimally emitting clouds. The failure of "simple" 
AGN photoionization models to descri be the narrow emissio n line 
spectrum of Seyfert 2 galaxies lead iBaldwin et alj d 19951) and 
iFerguson et al.l i 19971) to propose an alternative scenario in which 
the integrated narrow-line spectrum can be predicted by an integra- 
tion of an ensemble of clouds with a wide range of gas densities 
and distances from the ionizing source with an appropriate cover- 
ing factor and distribution. For each line, only a narrow range of 
density and distance from the continuum source results in maxi- 
mum reprocessing efficiency, corresponding to "locally optimally 
emitting clouds". 

• Dusty, radiation-pressure dominated photoionization models 
(Groves, Dopita & Sutherland |2004 Dopita et al. |2002|) . In these 
models, dust and the radiation pressure acting upon it provide the 
controlling factor in moderating the density, excitation and sur- 
face brightness of the photoionized gas. Additionally, photoelectric 
heat ing by the dust modifies the gaseous temperature structure. 

• iHumphrev etaD 12008) discuss in detail the need for a mixture 
of cloud properties (a range in U) in individual objects to explain 
the UV and optical rest-frame line ratios of high z radio galaxies. 



3.3 Stellar photoionization 

There are 3 objects, highlighted as blue triangles in Fig. 1, which 
show hybrid properties of both AGNs and HII galaxies. Their line 
ratios in diagrams a, b and c fit better a classification as HII galax- 
ies. Their lines are also relatively narrow compared with most ob- 
jects in the sample, for which the median value is 530 km s _1 while 
the hybrid objects have FWHM[OIII]<470 km s _1 . They show 
among the lowest [OII]/[OIII] and [OIII]/H/3 values in the sample. 

On the other hand, they show some features characteristic 
of active galaxies: they emit [NeV] and Hell (Fig. 1) and have 
large line luminosities (L[OIII]>10 42 erg s _1 ). The two objects 
for which [Nil] was measured occupy an intermediate region in the 
[OIII]/H/9 vs. [NII]/Ha diagram (Fig. 3) between HII galaxies and 
AGNs. We will call these hybrid objects. 

Based on the arguments exposed in the previous sections, 
AGN photoionization must play an important role in the ionization 
of the gas in type 2 quasars in this sample. The discrepancies of the 
standard AGN sequence show that a range of ionization and proba- 
bly physical properties must exist within the type 2 quasar sample. 
An internal range of cloud properties (e.g. density range, matter and 
bounded clouds, ionization level etc) is also likely to exist in indi- 
vidual objects. This is only natural, one cannot expect identical gas 
and continuum properties in all quasars, or ensembles of identical 
clouds in individual objects. As we discussed above, more sofisti- 
cated models that take this into account solve the problems of the 
standard AGN sequence. 

However, the emission line spectra of the three hybrid objects 
suggests that stellar photoionization might also be present with dif- 
ferent degrees of importance relative to AGN photoionization from 
object to object. Interestingly, although uncommonly high [OIII] 



luminosities for HII galaxies, the hybrid objects are at the lowest 
end of L[OIII] values within the type 2 quasar sample (<2x 10 42 
erg s _1 , while the median value is ~6x 10 42 erg s _1 ). This tenta- 
tively suggests that low line luminosities might be associated with 
relatively stronger stellar photoionization. 

We investigate next whether a contribution of stellar photoion- 
ized gas could solve the problems of the standard AGN models dis- 
cussed above. For simplicity, we have assumed that the total flux 
of a given line is due to the added contribution of the flux emit- 
ted by an AGN photoionized component (represented by the stan- 
dard AGN sequence) plus the flux emitted by a stellar photoionized 
component. To represent the spectrum emitted by t his gas, we have 
used t he real spectra of a variety of HII galaxies in lTerlevich et al.l 
Jl99ll) catalogue. We find from this study that the stellar ionized 
component must have [OII]/[OTII]>l in order to shift the models 
in the right direction in the diagrams. 

As an example we take UM448. The reason to choose this par- 
ticular HII galaxy is that it fulfills the [OII]/[OIII]>l requirement 
and has reported measurements of [OIII]A4363. Unfortunately, as 
for all other objects in the catalogue with [OII]/[OTII]>l, Hell is 
not measured. 

The flux of any line Flux tot relative to H/3* ot is given by: 



Flux to 
Hf3 tot 



Flux* + Flux 1 
H/3* + Hf3 AC 



where Flux* and H/3* are respectively the flux of the line 
and the flux of H/3 emitted by the stellar ionized gas. Flux AGN 
and H/3 AGN are the flux of the line and that of H/3 emitted by gas 
ionized by the AGN. 

Let us define x = h ^agn an d rearrange: 



Flux to 
H/3 tot 



+ X 



H/3" 



1 + X 



We assume [NeV]*/H/T << [NeV] AGAr /H/3 AGJV . When de- 
tected, HeII*/H/3* in HII galaxies is often in the range ~0.02-0.1 
(values as high as 0.4 are also possible, but are more typical of 
objects with very small [OII]/[OIII]<0.2 values). We will assume 
HeII*/H/3*=0.04. Using a different value within the expected range 
does not change our conclusions. 

We now create hybrid models by adding the stellar emission 
line spectrum to the AGN standard sequence of models and chang- 
ing the x value from sequence to sequence (see equation above). 
The ionization parameter U of the AGN photoionized gas (as in 
the standard AGN sequence) changes along each sequence. The re- 
sults are shown in Figs. 6, 7 and 8 as red dashed-lines. 

We find that by varying the relative contribution of the stellar 
to the AGN photoionized gas: 

• The discrepancy affecting the high and low ionization lines 
disappears (Fig. 6, panels e and /). 

• The large scatter observed in the HeII/H/3 diagram towards 
low values is now reproduced (Fig. 6, panel g). 

• The highest HeII/H/3 values (~0.8 for two objects) cannot be 
reproduced by the models (Fig. 6, panel g).To solve this problem 
with the stellar photoionized component, this should show similarly 
high HeII/H/3 ratios, which are extreme even for active galaxies. 
We notice this is also a problem for Groves, Dopita & Sutherland 
(2004) models, which require extremely high gas metallicities ~4 
Zq. No models in Ferguson et al. j 19971) reach these high HeII/H/3 
values either (see Fig. 3b in their paper). A strong contribution 
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Figure 6. As Fig. 1, adding the hybrid AGN+stellar models (red, dashed lines) for different x values shown in red (0.5, 1.0 and 4.0). The AGN sequence (black 
solid line) corresponds to x=0, i.e., the line fluxes are emitted by gas purely photoionized by the AGN. The arrows indicate the sense of increasing x. 



of matter bounded clouds to the integrated spectrum solves the 
problem (Binette et al. 1 19961) . Given the scarcity of these extreme 
HeII/H/3 ratios in type 2 AGNs in general, we wonder whether 
these two specific type 2 quasars are particularly rare. 

• Although the temperature problem is partially allevi- 
ated (Fig. 7), the hybrid models cannot explain the small 
[OIII]AA5007,4959/[OIII]A4363 measured in the quasar sample. 
Given the saturation of this line ratio at high electron temperatures 
(Osterbrock 1989), not even unrealistically high electron tempera- 
tures of a dominant stellar photoionized gas would solve the prob- 
lem. Notice that Binette et al. (1997) warned about the need for 
higher densities if the NLR dominates the emission, which would 
result in lower [OIII]AA5007,4959/[OIii]A4363 ratios for a given 
electron temperature. 

• The [NII]/Ha ratio is still too low compared with the obser- 



vations. Adding the stellar ionized component only makes things 
worse (Fig. 6, panel d and Fig. 8). However, as explained before, 
this could be suggestive of an overabundance of nitrogen. 

Another test for the stellar photoionization scenario consists 
of checking whether the measured continuum level is consistent 
with that expected for the stellar population responsible for H/3* . 
Let us consider the objects with the largest x values. This is the 
case of the hybrid objects, for which according to the exercise pre- 
sented here, x ~4 (see Fig. 6). For these three objects the Hf3 
luminosity is in the range ~3-6xl0 41 erg s _1 . We will assume 
L(H/?)=4,5 x 10 41 erg s~ 1 . If 80% is due to stellar photoionized gas 
then L(H/T )=3.6 x 10 41 erg s" 1 , rather similar to that of UMM448, 
L(H/3)=7.9xl0 40 erg s" 1 . The rest frame H/3 EW in the hybrid 
type 2 quasars is in the range ~32-38 A, implying EW values for 
H/3* ~25-30 A, quite similar to that measured in UM448 (43 A, 
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Figure 7. As Fig. 4, adding the hybrid AGN+stellar models (red, dashed 
lines) for different x values shown in red (0.5, 1.0 and 4.0). 
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Figure 8. As Fig. 3 with the hybrid models. 



Terlevich et al. Il99ll) . Therefore, for the hybrid objects the H/3* 
luminosities and equivalent widths expected from stellar photoion- 
ization are consistent with those measured for HII galaxies. 

We must not forget that the continuum emitted by type 2 
quasars is not necessarily stellar so that the measured level is an up- 
per limit to the stellar continuum level and the EW values discussed 
above for H/3* relative to the stellar continuum are lower limits. A 
more strict analysis would require a quantification of other possible 
contaminants, such as scattered light from the hidden AGN. Polari- 
metric information would be needed to characterize the nature of 
the continuum. Lacking this information, we can only say that the 
continuum level detected from the hybrid objects, where the stellar 
contribution is probably highest, is consistent with that expected 
from H/3*. 

While this is a simplistic exercise and one cannot expect all 
type 2 quasars to have identical spectra of the stellar ionized gas, it 
suggests nevertheless that adding a varying contribution of stellar 
photoionized gas works in the right direction to solve most of the 
problems affecting the standard AGN sequence. The temperature 
problem remains and a more sofisticated scenario with a range of 
gas densities or the presence of a matter bounded component might 
be a viable solution. 

Some studies suggest that the [OIIIJA5007 emission line is an 
unbiased indicator of the intrinsic optical-UV luminositiy of both 
type 1 quasars and radio galaxies. (Simpson 1998). According to 
the results above, it is possible that [OIII] has a strong stellar con- 
tribution in a fraction of type 2 quasars. Let us estimate the fraction 
of [OIII] flux originated by stars in the hybrid objects. We know 
that: 



[OIII] to 

H/3 tot 

[QUI]' 



[OIII] tot [OIII}* H/3* 
[OIII]* H/3* Hf3 tot 



using l ~^~pi i in the range 3.2-3.9 as measured for the hybrid 
objects, ^pl=2.9 (as forUM448) and -^^=0.8 (since x=4) and 
59-72%. This suggests that there could be 



H0* 

rearranging, { l OIII] L 



a fraction of type 2 quasars in which [OIII] is not a reliable indi- 
cator of AGN power. Other type of studies should be performed to 
investigate this issue more carefully (e.g. Simpson [l998l) . 



4 DISCUSSION AND CONCLUSIONS 

We have co mpared the optical em ission lines ratios of type 2 
quasars from Zakams ka et alj d2003h with standard AGN photoion- 
ization model predictions, Seyfert 2s, HII galaxies, and narrow line 
FRII radio galaxies. Moderate to high ionization narrow line radio 
galaxies and Seyfert 2s are indistinguishable from type 2 quasars 
based on their optical line ratios. The standard AGN photoioniza- 
tion models, valid for other type 2 AGNs, can reproduce success- 
fully the loci and trends of type 2 quasars in some of the main di- 
agnostic diagrams. On the other hand, these models are not exempt 
of problems and the discrepancies with the data are the same en- 
countered for other type 2 AGNs. The comparison between models 
and data suggests that a range of ionization and probably physical 
properties must exist within the type 2 quasar sample. An internal 
range of cloud properties (e.g. varying density) must also exist in 
individual objects. This is only natural, as one cannot expect iden- 
tical gas and continuum properties (e.g. ionizing luminosity) in all 
quasars, or ensembles of identical clouds in individual objects. 

Realistic models must take this into account. Possible solu- 
tions which have been extensively discussed in the literature are 
locally optimally emitting clouds, a mixture of matter and ioniza- 
tion bounded clouds, dusty, radiation-pressure dominated models 
or a mixture of clouds with different U values. 

The rele vant role played by AG N photoionization is not sur- 
prising, since Zakams ka et al.1 d2003l) selected objects with proper- 
ties characteristic of active galaxies. However, based on the lack 
of correlation between [OIII] and radio luminosities, other authors 
have suggested that an important fraction of type 2 quasars might 
not be do minated by AGN activity but by star formation (Vir Lai & 
Ho l2007l : see also Kim et al. l2006h . 

We have found that stellar photoionization is obvious in a 
small fraction of objects (3 out of 50) which are characterized by 
low [OIII] luminosities and large [OII]/[OIII] ratios compared with 
mos t type 2 quasars in the s ample. 

IZakamska et all d2003h sample is strongly biased towards ob- 
jects with high line luminosities ( L[OIII]>10 42 erg s" 1 ). L[OIII] 
can be as low as ~10 40 erg s _1 in radio-quiet type 1 quasars and 
na rrow l ine radio galaxies (e.g. Bennert et al. 2002, Tadhunter et 
al. 1998). There must be many type 2 quasars with L[OIII] in the 
range ~10 40-42 erg s _1 . The hybrid objects discussed here, where 
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stellar photoionization is obvious, have a mong the lowest L[OIII] 
values within the type 2 quasar sample of IZakamska et al.l J2003h . 
This tentatively suggests that stellar photoionization could be rela- 
tively more important in type 2 quasars with lower [OIII] luminosi- 
ties. If these objects have lower power AGNs (e.g. Simpson l 19981) . 
it would seem plausible to find a relatively higher contribution of 
the stellar photoionized gas to the observed emission line spectrum. 
The fraction of hybrid objects, therefore, could be much larger, if 
lower L[OIII] values were considered. 

Since star formation has been found in differents classes of 
type 2 AGNs at different z (e.g. Holt et al. 120071 Tadhunter et al. 
1 20051 Alonso-Herrero et al. 2008j) a nd in at least a fraction of type 
2 quasars (Martfnez-Sansigre et al. 120081 . Lacy et al. 2007), stel- 
lar+AGN photoionization is a plausible scenario. Inspired by these 
arguments, we have explored an alternative scenario to pure AGN 
photoionization in which a varying contribution of stellar ionized 
gas is added to the line fluxes in type 2 quasars. Although the hybrid 
models presented here are rather simplistic and no reliable quanti- 
tative results can be extracted about the relative importance of stel- 
lar vs. AGN photoionization, they reproduce the type 2 quasar ra- 
tios (and type 2 AGN in general) quite successfully (better than the 
standard AGN sequence). This suggests that stellar photoionization 
might also be present in many type 2 quasars, in addition to AGN 
photoionization. Given that other type 2 AGNs have similar line 
ratios, this applies as well in those cases. On the other hand, and 
contrary to the more sofisticated AGN models discussed in the lit- 
erature, the hybrid models cannot solve the "temperature problem" 
(see §3.1 B). Regarding all other line ratios, given the strong de- 
generacy with more sofisticated AGN photoionization models, it is 
not possible to favour one scenario or another in terms of the line 
ratios. 

Other sources of information would be very valuable to test 
whether stellar photoionization is present in type 2 quasars and 
characterize its importance relative to the AGN photoionization 
(e.g. evidence for extended star formation, possible correlations be- 
tween star formation indicators and the line ratios which require a 
stronger contribution of stellar ionized gas, lower continuum po- 
larization level in objects with hints of stellar photoionization from 
the line ratios, etc). If the emission lines, in particular [OIII]A5007 
have a strong contribution of stellar photoionized gas, the [OIII] 
line might not be a good indicator of AGN power. 
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